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The utilization of satellites for communication systems has expanded 
considerably in recent years. C and Ku-bands of frequencies are already 
congested because of high demand. Future directions of satellite 
communications are moving towards Ka and V-bands. Earth to satellite 
communications are moving towards higher frequency bands in future which 
are more sensitive to environment. Rain causes severe degradation in 
performances at higher frequency bands specially in tropical regions. Several 
mitigation techniques are proposed to design reliable system. Time diversity 
is one of the potential candidate for it. However, time diversity analysis 
requires measured rain attenuation data. For future high frequency link 
design those data are not available at most of the places. This thesis proposes 
a method to utilize 1-minute rain rate to analyze time diversity technique at 
any desired frequency. This paper proposes a method to utilize 1-minute rain 
rate to analyse time diversity rain rate gain. In proposed method, it is 
assumed that rain rate gain with delay can represent rain attenuation gain 
with delay for same period of time at same location. The characteristics of 
rain rate and rain attenuation almost same because the attenuation causes due 
to rain. One year measured rain rate in Malaysia is used to predict rain rate 
gain. The measured gain at 12.225 GHz signal is compared with that 
predicted by ITU-R based on rain rate measurement and is found good 
agreement. Hence it is recommended that the time diversity gain can be 
predicted using measured rain rate for any desired frequencies. 
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1. INTRODUCTION 


Earth-to-satellite links that operate at frequencies higher than 10 GHz (Ku, Ka, and V bands) are 
strongly affected by propagation impairments especially by precipitation the tropical climate and equatorial 
climate experience heavy precipitation all over the year [1]. Therefore, the satellite communication system 
performance suffers from severe degradation at high frequencies [2]. It is not practical to include very high 
fixed fade margins in link budgets as a compensation for the great propagation impairments due to effects of 
transmitter power limitations. In order to overcome such problems, implementation of Fade Mitigation 
Techniques (FMTs) becomes significant to ensure the continuity of system availability [3]. 

FMTs have to take in consideration the system performance aims, the system architecture 
parameters and the frequency bands in operation. Time diversity is one of the best choice among FMTs in an 
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Earth-satellite links to encounter severe channel fading in terms of cost and benefits [4]. The performance of 
time diversity is generally evaluated using measured rain attenuation time series with and without delay of 
time. It calculates the joint exceedance of rain attenuation probability both with and without delay [5]. 
Time diversity gain is the indicator for improvement. An empirical model to predict time diversity gain 
function was derived in Malaysia based on measured data to estimate rain attenuation exceedance 
probability [6]. Similar gain model as a function of delay time, attenuation level and frequency in operation 
was proposed in [7, 8], two methods were developed for calculating joint exceedance rain attenuation 
probability which represent the complementary cumulative distribution function of rain attenuation that 
follows the lognormal distribution 

In rain attenuation analysis along earth-to-satellite links, it can be observed obviously that 
attenuation of signal is directly related with intensity of rain occurred [9]. Hence the variations of 
attenuations must follow the variations of rain intensity along the propagation path [10-12]. onsidering this 
fact in mind, real time rain intensity variations can be utilized to evaluate time diversity improvement on 
earth-to-satellite links which is conventionally estimated using real-time rain attenuation data [13]. 
This paper investigates the possibility to use rain rate time series to estimate time diversity gains for earth-to- 
satellite links [14-17]. Measured rain rate and rain attenuation data at Ku-band are used for investigation. 

The equation of the marginal distribution which is converted conventional cumulative probability 
distribution function (CDF) is shown in below [7] 


P(A) = P[A(t) > A] = f° E[A()]dA (1) 


where €(A) is the density function of rain attenuation A(t), and P(A) is its integral, calculated for all time t of 
the rain attenuation time series. The joint distribution is given by [7] 


P[A(t) > A,ACE+T) >Al= SP SP y[AW, A(t + T)]dA(Od A(t + T) (2) 


where y[A(t), A(t + T)] is the joint probability density function of A(t) and A (t +T) with delay time T. Most 
of Time Diversity studies follow | to 60 minutes delay of times as shown in Figure I(a) and 1(b). 

The complementary cumulative distribution function of rain attenuation conditioned to the time 
delay can be used to estimate the performance of time diversity as follows: A(t) refers to rain attenuation 
without time delay, while T is defined as a time diversity delay and A (t +T) refers to rain attenuation 
shifted of T. 

This represents the time of percentage exceeded for minimal rain attenuation and the signal of 
satellite to earth will be subjected to minimal rain attenuation as a result of using time diversity technique and 
is expressed by the model 


Arp = min[A(t), A(t + T)] (3) 


There is a need to make a comparison of the original time series of attenuation and the same time 
series delay by given time interval in order to get the lowest attenuation value as explained in (3). 
An illustrative example for the use of time delay technique is presented in Figure 1. An event of measured 
signal on the satellite-Earth link is applied. It is clear that the maximum attenuation of the signal without time 
delay is about 28 dB, but after applying time diversity technique with time delay of 15 and 30 minutes, the 
maximum rain attenuation is 24 and 7 dB respectively. 

The signal delays T-minute or in the range 1-60 minutes instantaneously. From delaying the signal 
may be written as the probability from the diversity combination experiences an attenuation 
greater than A [11]. 


prob(Agiy = A) = 1-PI(A, or Ar) < A] (4) 
The time-diversity gain (in dB) for a specified exceedance probability (p) as A,/ (non-diversity). 
The cumulative distributions of rain attenuation A, with different time delays were analysed based on 


measurements done in Malaysia [12] are shown in Figure 2. From figure, it is clear that more than 20 dB 
improvement is possible through time diversity technique. 
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Figure 1. An Example, which illustrates the use of time diversity technique and its improvement in the signal 
of a satellite-Earth link for an event, measured at Penang-Malaysia on October 21, 2009, 
(a) Delay time T=15 minutes, (b) Delay time T=30 minutes [6] 
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Figure 2. Cumulative distributions of rain attenuation with time diversity for delays in the range of 1-60 
minutes measured at 12 GHz in Malaysia [12] 


Time diversity gain is the value of the difference in dB between rain attenuation of the 
communication signal and the same signal with defined time delay [18]. Therefore, the time diversity gain is 
used to evaluate the time diversity improvement in performance as a fade mitigation technique for rain 
attenuation [14]. Time diversity gain can be expressed instantaneously by the following equation. 


GArp = A(t) — Arp (dB) (5) 


where GA 7p is time diversity gain using rain attenuation. 

The difference at equal probability, time diversity statistical gain can be described by the rain 
attenuation A;yp) exceeded from the marginal distribution [19] as the long-term annual probability 
distribution and A(7),«,p) of the joint distribution relative to delay T: 


GArpmp) = Ayp — A(T ogy (6) 


However, time diversity analysis requires measured rain attenuation data. For future high frequency 
link design those data are not available at most of the places [20-21]. This paper proposes a model to utilize 
1-minute rain rate to predict time diversity gain at any desired frequency. In proposed method, it is assumed 
that rain rate with delay can represent rain attenuation with delay for same period of time at same location. 
This assumption is valid as long as the attenuation causes due to rain. One year measured 1-minute 
integration rain rate at International Islamic University Malaysia is used for analysis. 
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2. RAIN RATE DATA COLLECTION 

The rain rate data were collected for 1-year 2014, from January | to December 31 at International 
Islamic University Malaysia (IUM) Kuala Lumpur campus (3.257° N, 101.7325 °E). A real-time rain gauge 
of Casella type with 0.2 mm bucket size was used to measure the data. The data collection setting was 
10 second integration time. Readings were displayed by a personal computer, which was programmed to 
convert the data in 1- minute integration time. 


3. PROPOSED TIME DIVERSITY GAIN USING RAIN RATE 

Rain intensity and attenuation caused due to rain are closely related to each other. Based on the rain 
attenuation equation with time delay, the complementary cumulative distribution function of rain rate with 
time delay is proposed to be used to estimate the performance of time diversity as follows: 

R(t) refers to rain rate without time delay, R (t+7) refers to rain rate shifted of T. while T is defined 
as a time diversity delay. (7) to (9) are proposed using rain rate time series as follows: 


P(R) = P[R(t) > R] = f, &[R(t)]dR (7) 


where €(R) is the density function of rain rate R(t), and P(R) is its integral, calculated for all time t of the 
rain rate time series. The joint distribution is proposed as 


P[R(Q)>R,R(C+T)>RI= fe fe y[RWO,R(t + T)]dR(e)dR(t + T) (8) 


where y[R(t), R(t + T)] is the joint probability density function of R(t) and R (t+T). Hence the rain rate with 
time diversity can be expressed as follows 


Rrp = min[R(t), R(t + T)] (9) 


Conventionally used time diversity gain is the value of the difference in dB between rain attenuation 
of the communication signal and the same signal with defined time delay. The proposed time diversity gain 
using rain rate can also be defined instantaneously as the difference of rain rate with and without delay in 
mmn/hr using (6) as follows: 


GRrp=R(t)-Rrp(mm/hr) (10) 


The difference at equal probability, time diversity statistical gain can be found by the rain rate R,,.~) 
exceeded in the marginal distribution as the long-term annual probability distribution and Ry, of the joint 
distribution relative to delay T: 


GRrpcmp) = Rmxp — Roop (11) 


This Figure 3 shows the rain rate from different percentage which was derived by using the program 
and was transmitted twice with time delay, T. This Figure 3 shows that rain rate decreases as T increases. 
Thus, the time delay with several minutes especially in small cumulative time percentage is shown large 
diversity gain in the Figure 4. Most probably this method is powerful in future for rain fade mitigation. 
In future this technique will be demandable for the transmitter and receiver. This technique can be used by 
lower cost in the system for millimetre band satellite systems. If time diversity with time delay, T, effective 
rain rate can be reduced 10, 15, 30, 40, 52 and 68 mm/hr at 1, 3, 5, 10, 20 and 30 minutes delay even at 
cumulative time percentage at 0.01%. Variation in cumulative time percentages as a function of delay 
time, T. 
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Figure 3. Proposed complementary cumulative distributions of 1-minute rain rate measured at IIUM with 
several time delays for the year 2014 
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Figure 4. Variations of proposed time diversity rain rate gain in mm/hr for different percentages of time with 
several time delays 


4. PROPOSED TIME DIVERSITY GAIN USING PREDICTED ATTENUATION 

Based on measured rain rate at 0.01% without and with time delay are estimated from Figure 4 for 
0, 1. 3, 5, 10, 20, 30 minutes and found 125, 115, 108, 97, 85, 73, 66 mm/hr respectively. The cumulative 
distributions of rain attenuation A by considering Roo1q% as 125, 115, 108, 97, 85, 73, 66 mm/hr rain rates are 
predicted by ITU-R equations (1 to 9) [17] and shown in Figure 5. In prediction, International Islamic 
University Malaysia Kuala Lumpur campus is considered as receiver location and MEASET-3 as satellite 
with 12.255 GHz frequency for transmission and an elevation angle of 77.4°. From figure, it is clear that 
more than 12 dB improvement are possible through time diversity technique. For illustration, at 0.01%, the 
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rain attenuation level without applying time diversity technique is 29.4 dB. While, the rain attenuation levels 
are 27, 25.5, 24, 21.5, 19.6, 15.6 dB diversity with time delays of 1, 3, 5, 10, 20 and 30 min respectively. 
Figure 6 shows strong similarities with results. 

Based on (6), the time diversity gain is estimated using Figure 6 with 0.1%, 0.01% and 0.001% and 
shown in Figure 7. Table | presents the gain with 1, 3, 5, 10, 20 and 30 minutes delays. It is obvious that gain 
increases with longer delays for all percentage of time. For same period of delay, gain at lower outages are 
always higher than that of higher percentages. These findings are found very similar with those obtained by 
measurements and shown in Figures 5 and 6. 
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Figure 5. Complementary cumulative distributions of predicted rain attenuation for satellite-earth link at 
12.225 GHz with 0, 1, 3, 5, 10, 20 and 30 minutes time delays 
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Figure 6. Predicted time diversity rain attenuation gain for different percentages of time with several 
time delays 
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The time diversity gain was estimated using measured rain attenuation distribution from the 
MEASET-1 satellite signal transmitted at 12 GHz at the KualaLumpur, Malysia with an elevation angle of 
77.4° [18]. The measured gain at 0.01% is compared with those predicted in Table | with delays of time 1, 3, 
5, 10, 20 and 30 minutes and presented in Figure 7. It is obvious that the gain is almost similar until 10 
minutes delay and predicted gain is 1.6 dB higher for 30 minutes delay. 
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Figure 7. Comparison Measured time diversity rain rate gain and predicted rain attenuation gain for different 
percentages of time with several time delays 


Table 1. Predicted Rain Attenuation Gain at 12.225 GHz 








Top 

Time 0.1 0.01 0.001 
1 1 18 25 
3 1.5 2.9 5 
5 2.3 3.8 6 
10 3.5 6 8.5 
20 5 9.7 13.5 
30 6.3 12.4 17.5 





Hence it is recommended that the time diversity gain can be predicted using measured rain rate with 
very good accuracies and its shown in Figure 8. Since, measured |-minute rain rate is available in almost all 
locations of the world, it can be utilised to predict time diversity gain for any desirable frequencies. 
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Figure 8. Comparison between the measured attenuation gain with those predicted at time delays of 1, 3, 5, 10, 20 
and 30 minutes for 12 GHz 


5. CONCLUSION 

Time diversity is one of the techniques proposed for mitigation of rain fade. It is a clear, the technique 
of time diversity is suitable for applications such as data transfer process that not needed until the rain fades 
over. However, time diversity analysis requires measured rain attenuation data which are not available at 
most of the places. This paper proposes a method to utilize 1-minute rain rate to analyse rain rate gain which 
is used for time diversity gain at any desired frequency. It is assumed that rain rate gain with delay can 
represent rain attenuation gain with delay for same period of time at same location. Measured 1-year rain rate 
distribution with and without delay were used for rain rate gain. Measured rain rate is used to predict rain 
attenuation distributions by ITU-R method. This distribution with and without delay is found very similar 
characteristics with those obtained by measurements. Predicted rain attenuation distribution is analysed to 
estimate attenuation gain and is compared with measured attenuation gain for 12.225 GHz. The predicted 
attenuation gain is found close to measured one. Since, measured 1-minute rain rate is available in almost all 
locations of the world, it is recommended to be utilised to predict time diversity gain for any 
desirable frequencies. 
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